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ABSTRACT 

We present R ~ 600, lO/mi - 37/im spectra of 53 Ultraluminous Infrared 
Galaxies (ULIRGs) at z < 0.32, taken using the Infrared Spectrograph on board 
the Spitzer space telescope. All of the spectra show various fine structure emis- 
sion lines of Neon, Oxygen, Sulfur, Silicon and Argon, as well as one or more 
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molecular Hydrogen lines. Some objects also show emission lines from Chlorine, 
Iron, Phosphorous and atomic Hydrogen, as well as absorption features arising 
from C 2 H 2 , HCN, and OH . We employ diagnostics based on the fine-structure 
emission lines, as well as the equivalent widths and luminosities of polycyclic 
aromatic hydrocarbon (PAH) features and the strength of the 9.7/im silicate ab- 
sorption feature (S s u), to explore the power source behind the infrared emission 
in ULIRGs. We show that the infrared emission from the majority of ULIRGs is 
powered mostly by star formation, with only ~20% of ULIRGs hosting an AGN 
with a comparable or greater IR luminosity than the starburst. The detection of 
the [NeV]A14.32 line in just under half the sample however implies that an AGN 
makes a significant contribution to the mid-IR flux in ~ 42% of ULIRGs. The fine 
structure line ratios, luminosities and PAH EWs of our sample are consistent with 
the starbursts and AGN in ULIRGs being more extincted, and for the starbursts 
more compact, versions of those in lower luminosity systems. The excitations 
and electron densities in the narrow-line regions of ULIRGs appear comparable 
to those of lower luminosity (10 10 <Lj r (L Q ) < 10 115 ) starbursts, though there is 
evidence that the NLR gas in ULIRGs is more dense. We show that the combined 
luminosity of the [NeII]A12.81 and [NeIII]A15.56 lines correlates with with both 
infrared luminosity and luminosity of the Q.2fim and 11.2/mi PAH features in 
ULIRGs, and use this to derive a calibration between PAH luminosity and star 
formation rate. Finally, we show that those ULIRGs with 0.8 < S s u < 2.4 are 
likely to be powered mainly by star formation, but that those with S si i < 0.8, 
and possibly those with S s u > 2.4, contain an IR-luminous AGN. 

Subject headings: infrared: galaxies — galaxies: active — galaxies: starburst - 
galaxies: evolution 



1. Introduction 



Ultraluminous Infrared Galaxies (ULIRGs, those objects with 1-1000 /im luminosities in 
excess of 1O 12 L ) were first discovered in the 1970s (IRieke fc Lowlll972l ). Since then, they 
have fascinated astronomers with their unique and extraordinary properties, and infuriated 
them with their singularly opaque natures, almost in equal measure. 

In the local Universe ULIRGs are a rare, if interesting oddity, with only fifty or so 
examples known at z < 0.1. First un covered in significant numbers by su rveys with the 
Infrared Astronomical Satellite (IRAS, ISoifer et all Il984l ; iHouck et al.lll985l ). work focused 
on determining the power source behind their colossal infrared emission. This initially pro- 
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voked heated debate between a 'starburst' camp and an 'AGN' camp, a debate that has 
not yet entirely cooled. Early s tudies showed that som e ULIRG optical spectra resem- 
bled those of starburst galaxies Jjoseph fc Wright 1985h . whereas others contained emis- 
sion lines characteristic of Seyferts ( Sanders et all 1988H . Radio observations showe d di- 
rect evidence fo r starbursts in some ULIRGs (ICondon et all Il99ll ; ISmith et all Il998l ) and 
AGN in others (ILonsdale et al.l 120031 ; iNagar et all 120031 ) . Over the last decade or so how- 
ever a consensus has started to emerge; local ULIRGs are likely to be 'composite' objects, 
with most powered mainly by a starburst, but with a significant fraction also containing 
an IR-luminous AG N. This is suggested from several lines of eviden c e, including opti- 



cal/UV spectrosco py (jVeilleux et al 



1995; Vei 



1998; 



Farrah et al.ll2005l). mid-infared spectroscopy ( 



Rigopoulou et allll999l ; iGenzel &: Cesar skv 200 01; iTran et al.ll200lD. mode lling of their 



l-lOOO/an sp ectral energy distributions dKlaas et al.l l200ll ; iFarrah et al.l 120031 ) . and X-ray 



leux. Sanders & Kim 


1999; I 


jipari et al. 


2003; 


Lutz et al. 


1996; 


Genzel et a 


1998; 


Lutz et al. 



observations ( Franceschini et al 2003 ; Ptak et al 20031). though some r ecent studies suggest 



y wi th galaxy mergers (IFarrah et al 



a significantly greater average AGN contribution (llmanishi et al.l 12007). Local ULIRGs are 
also associated almost exclusive 



2001; Bushouse et al 



(Sanders et al 


1988: 


Tacconi et al 


2002; 



2002 ), and may be involved to some degree in triggerin g 



Excellent reviews of the properties of ULI RGs can be found in 
and more recently in Lonsdale et al. (j2006 ). 



Zauderer et a 



Sanders fc Mirabel 



2QQ7J). 



(1996) 



Their rarity in the local Universe compared to lower luminosity systems initially led 
astronomers to believe that ULIRGs did not play a fundamental role in galaxy formation 
processes, however this perception changed abruptly when it was realised that ULIRGs 
were vastly more numerous at high r e dshift. First hinted at by spectroscopic followup 



of IRAS surveys ([Hacking et al.l Il987l ; ILonsdale et al.l Il990l ; ISaunders et all 119901 ) , which 
showed strong evolution in the ULIRG luminosity function with reds hift, and from th e 



discovery of a remarkably high cosmic inf rared background by COBE ( Puget et al 
this was confirm e d by surveys with ISO (IRowan-Robinson et al.l Il997l ; iLevine et al. 



996 



1998 



Dole et al.l l200ll ; IVerma et al.l l2005l ). which found a large population of ULIRGs up to 



z ~ 1.5, and thrown into sharp relief by sub- mm surveys (IHughes et all Il998l ; lEales et al 



2000l ; iBorys et all 120031 ; ICoppin et al.l 120061 ) , which showed that there were several hundred 
ULIRGs per square degree at z > 1. Though obviously much harder to study, these dis- 
tant ULIRGs seem superficially similar to their low redsh ift counterparts in that they ap 
pear to be powered by both starburst and AGN activity (IFarrah et al.l l2002bl ; ISmail et al. 



2003, 2004; lAlexander et al. 



2005 



Takata et al 



2006 



Valiante et al.l 120071 ). and are proba- 



bly mergers (IFarrah et al.ll2002al ; IChapman et al.ll2003l ). Their properties may make them 
important tools in understanding the global evolution of galaxies and large-scale structures; 
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their rapid star formation rates and comoving number densities make them strong can- 



didates for being the rapid growth phases of massi ve elliptical galaxies (jScott et al.l 12002 



ap 

l. : 



Rocca-Volmerange et al.ll2004l ; ISwinbank et al.ll2006f). and they may serve as efficient 'light - 



houses' of the seeds of massive clusters at z > 1.5 (jBlain et al.ll2004l ; iFarrah et al.ll2006a| jbl). 



The still controversial nature of the power source in local ULIRGs, coupled with the 
central position that ULIRGs seem to play in several astrophysical processes at high redshift, 
makes it ever more important to understand the nature of the heavily obscured starburst and 
AGN activity in these systems. This is ideally done in the mid-infrared, directly sampling 
the emission from the hot dust that shrouds the central engines of ULIRGs, and using 
fine-structure lines, which suffer much less from exti nction effects than optical or near-IR 
lines. The recently launched Spitzer space telescope (jWerner et al.l 120041 ) provides an ideal 
platform to undertake such studies, wi th its suite of mid-IR instruments, including the 
Infrared Spectrograph (jHouck et alll2004j ). which offers dramatic improvements in sensitivity 
and resolution compared to previous generation facilities. In this paper, we present high 
resolution mid-infrared spectra of 53 local ULIRGs, and discuss some spectral diagnostics 
based on their emission line fluxes and other spectral features. We assume a spatially flat 
cosmology with H = 70 km s _1 Mpc -1 , ft = 1, and fl m = 0.3. 



2. Analysis 
2.1. Observations 

The ULIRGs presented here were observed as part of a large study within the IRS 
GTO program to obtain mid-infrared spectra of 110 low redshift ULIRGs (Spitzer pro - 



gram ID 105). These 110 UL IRGs were selected from the IRAS Uy ( ]Kim fc Sanders! Il998l ) 
and 2Jy (IStrauss et al.l Il990l ) spectroscopic surveys, and from the FIRST ULIRG sample 



( jStanford et al.ll2000l ). Low resolution spectra were obtained of all 110 objects, and high res- 
olution spectra were obtained of the 53 brightest (at 60/iin, those with / 60 > 0.7Jy) objects. 
The low resolution spec tra span 5.2iim - 38 .5 /xm, with a r e soluti on of R~ 60 — 12 5. Initial 



results are presented in lArmus et al.l (120041 ); ISpoon et al.l (20041) and lArmus et al.l (120071 ). 
Molecular Hydrogen masses are presented in Higdon et al.l (120061 ). and crystalline silicate 
measurements are presented in ISpoon et all hood ). An atlas of the low resolution spec- 
tra can be found in iDesai et al.l (120071 ). including Polycyclic Aromatic Hydrocarbon (PAH) 
luminosities and equivalent widths. Measurements of the strengt hs of the 9.7 u m sil icate 
absorption features from the low resolution spectra can be found in ISpoon et al.l (120071 ) . 



Here, we present the 53 high resolution spectra. The sample is listed in Table HJ A few 
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of the sample have IR luminosities that lie slightly below the canonical ULIRG lower limit 
of 10 12 L Q , but for simplicity we refer to them as ULIRGs for the remainder of this paper. 
Each ULIRG was observed with both the Short-High (SH, 9.9/mi - 19.6/im, 11.3" x 4.7", 
R~ 600, 2.3" pix- 1 ) and Long-High (LH, 18.7/^m - 37.2/mi, 22.3" x 11.1", R~ 600, 4.5" pix" 1 ) 
modules onboard the IRS. The targets were placed in the center of each slit by performing 
'high' accuracy peak-ups using the blue peak-up array, on either a nearby 2MASS star or on 
the nuclei of the ULIRGs themselves, and observed in two nod positions. For five ULIRGs, 
their nuclei are separated by distances of ~5" or more; for IRAS 08572+3915 the slits were 
centered on the north-western nucleus, for IRAS 14348-1447 the slits were centered on the 
south-western nucleus, for IRAS 19254-7245 the slits were centered on the southern nucleus, 
for IRAS 23498+2423 the slits were centered on the north-western nucleus, and for Mrk 463 
the slits were centered on the eastern nucleus. The other ULIRGs are all either single-nucleus 
(in ground-based imaging) systems, or have close separation (<5") double nuclei. While the 
available optical spectroscopy for our sample is not homogenous, there are no examples of 
objects in our sample where two nuclei with clearly different optical spectral classifications 
fall within the IRS slits. 

For most of the SH observations we observed each object for 6 ramps, with a ramp time 
of 30s, to give a total on-source exposure time of 180s. For the fainter targets we observed 
for two ramps using a ramp time of 120s, for a total on-source exposure time of 240s. The 
LH observations were the same for each object; namely 4 ramps and a ramp time of 60s, for 
a total on-source exposure time of 240s3. 



2.2. Data reduction 



The data were processed through the Spitzer Science Center's pipeline reduction soft- 
ware (version 13.2), which performs standard reduction tasks such as ramp fitting and dark 
current subtraction. To ensure an accurate flatfielding correction we started our reduction 
from the unflatfielded (droopres) images. Starting with these frames, we flagged rogue and 
otherwise 'bad' pixels using the 2rsc/ear|^| tool, which uses a mask of rogue pixels for each 
campaign to first flag and then replace rogue pixels. The individual frames were then com- 
bined into a singl e image, and spect ra were extracted from each nod using the SMART 
software package (IHigdon et al.l 120041 ). using full-slit extraction. Wavelength and flux cali- 



1 Details for each observation can be found by referencing the AOR keys given in Table [T] within the 
Leopard software, available from the Spitzer Science Center 

2 This tool is available from the SSC website: |http: / / ssc.spitzer.caltech.edu 
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bration were performed by dividi ng the extracted sp ectra by that of a standard star, £ Dra, 



and multiplying by its template (jCohen et al.l 120031 ) . Features present in only one nod po- 



sition were treated as artifacts and removed manually. The two nods were then combined. 
The pixels on the edge of each order (typically the first and last 12 pixels) corresponding to 
regions of decreased sensitivity on the array were then removed to give the final spectrum for 
each object. The resulting spectra were generally of excellent quality. In a few cases some 
slight mismatch in continuum fluxes between orders was apparent, but not to an extent that 
could impact the analysis of emission line fluxes. 

Both the SH and LH slits are too small to allow for on-slit background subtraction. 
Ideally the sky continuum background (which is comprised mainly of temporally and spatially 
varying zodiacal light) should be subtracted using contemporaneous 'sky' observations, taken 
with similar exposure times and as close on the sky as possible to the observations of the 
target. Such observations were however not ta ken for our sampl e. In the absence of sky 



observations, modelled sky fluxes can be used (jReach et al.ll2003l ). but they are uncertain 
by at least a factor of two. In this paper we are interested only in the emission line fluxes, 
for which background continuum subtraction is not necessary. Hence, we do not correct our 
spectra for contamination from sky continuum background. 

We measured line fluxes and wavelengths by assuming the emitting region was a point 
source at the spatial resolution of the SH and LH slits, subtracting the continuum via 
a cubic spline fit over a ~ 0.5/im region centered on each line, and then fitting a single 
Gaussian profile to each line. We found a pure Gaussian profile to be a good fit to the lines 
in virtually all cases. The resulting 3<r uncertainties in the wavelengths are themselves a 
function of wavelength, and are typically 0.01/xm at 10/zm, 0.03/mi at 20/im, and 0.04/zm 
at 30/im. Blended lines were measured by simultaneously fitting multiple gaussians to the 
combined profile. Upper limits were determined by measuring the noise level of the data at 
the wavelength where the line is expected to lie; this can give rise to significant variations in 
upper limits for different lines in the same source if the wavelength of a line lies in an order 
overlap region. 



2.3. Extinction corrections 



Previous authors have generally used either near-IR line ratios, or the [SIII]AA18. 713, 33.481 



line ratio to derive an extinction correction (e.g. iVerma et al.ll2003l ). however we only detect 
both [SIII] lines in a small fraction of our sample, and the available near-IR spectroscopy 
for our sample is heterogenous and incomplete. Given this, and the uncertain structure of 
the narrow-line emitting gas in the mid-IR, we have chosen not to correct our line fluxes for 



- 7- 



extinction. We can however estimate what effect this lack of an extinction correction will 
have on our analysis. The relationship between the intrinsic and observed flux ratio for a 
pair of emission lines at wavelengths Ai and A2 can be written: 



^2 / int 



I_M 



x 10 



0.4 



obs 



where A\ 1 , A\ 2 and Ay are the extinctions at the wavelengths of the pair of lines and in 
the rest-frame V band, respectively. Assuming Ay = 3 .169 x E(B — V) , and obtaining 
values for A\/E(B — V) from a standard extinction law (ILi &: Draineil200ll ) then allows us 
to estimate the effect on a line ratio in terms of an increase in the V band extinction. Using 
this formalism, we find that the effect on most line ratios if Ay is increased is small. For 
example, for the [NeIII]/[NeII] ratio we obtain: 



V \NeII\ j M \ [Nell] J ^ 1 > 

which is negligible unless the increase in Ay is at least a few tens. Hence our lack of an 
extinction correction, while important to be aware of, should not unduly affect our analysis 
of most line ratios, though the effect of a lack of an extinction correction on line luminosities 
is more significant. We discuss the effects of extinction in terms of an increase in Ay for 
both fine structure line ratios and luminosities in the following sections. 

Finally, we note that the magnitude, and in some cases the direction, of the effect on 
a line ratio for an increase in extinction depend significantly on ones choice of extinction 
law. In Table [2] we present the scaling factors for several line ratios for an increase in 
Ay, assuming some commonly used extinction laws. For a given line ratio there is broad 
consistency between the extinction laws, but differences of up to 25% between the scaling 
factors are common. These differences in scaling factors between extinction laws should be 
kept in mind in the following sections. 



3. Results 



The SH spectra are presented in Figures [H to El and the LH spectra are presented in 
Figures El to EH 
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3.1. Common lines 



All of the spectra show various fine structure emission lines of Neon, Oxygen, Sulfur, 
Silicon and (depending on redshift) Argon. Also present are molecular Hydrogen lines. 
These lines are listed in Table |3j We detect [NeII]A12.81 and [NeIII]A15.56 in nearly every 
object. [SIII] A18.71 is also common, detected in ~80% of the sample, and [SIV]A10.51 is 
detected in just under half of the sample. The detection of [SIII] A33.48 is dependent on 
redshift, requiring z < 0.11 (for the line to lie at < 37/im in the observed frame), but 
is detected in ~80% of the objects where this line lies within the LH bandpass. Other 
lines whose detection is dependent on redshift include [Aril] A8. 99, which is detected in the 
majority of objects at z > 0.14, and [Sill] A34.82, which is seen in about half of the objects at 
z < 0.06 (though this line lies in a noisy LH order even at z = and thus is hard to detect). 
Three higher ionization lines are also present, though less commonly than the lines discussed 
above; [NeV]A14.32 and/or [OIV]A25.89 are detected in just under half the sample, while 
[NeV]A24.32 is detected in about one third of the sample. Turning to molecular Hydrogen 
lines; the S(3), S(2) and S(l) pure rotational transitions of H 2 are seen in nearly all of the 
sample, while the S(0) H 2 transition is seen in about one third of the sample, though its 
rarity compared to the other molecular Hydrogen lines is probably as much to do with the 
rising continuum toward s longer wave l engths as an ythin g else. Our line fluxe s are in all cases 
consistent with those in lArmus et al.l (120041 . 120071 ) and iHigdon et al.l (120061 ). though we use 
a more recent version of the IRS pipeline. The fluxes reported here should therefore be more 
accurate than those previously published. 

In addition to fine-structure (and other) emission features, some objects show one or 
more absorption features. The focus of this paper is the fine-structure emission lines, so 
here we only briefly mention these features, deferring a full discussion to later work. Several 
objects show two absorption features, corresponding to the vibration-rotation absorption 
bands of C 2 if 2 a t 13.70/xm an d HCN at 14.02/im. An extensive discussion of these features 
can be found in lLahuis et al.l (120071 ). A total of ten objects, most prominently Arp 220, 
IRAS 15250+3609 and IRAS 20551-4250, show an OH absorption feature at ~34.6/im (IRAS 
15250+3609 may also show a further OH absorption feature at rest-frame 28.9/xm, though the 
significance of detection is weak). This feature is most likely the 2 n 3 / 2 J = 3/2— 2 Hi/ 2 J = 5/2 
OH absorption doub let, which is thought to pump the 1667MHz OH maser line, among others 



felitzur et aDll976h. This absorption f eature has been seen p reviously in G alactic sources 



(e.g. lJusttanont et al.lll996f ). NGC 253 (IGoicoechea et al.ll2005h . and Arp 220 (ISkinner et al. 



19971 ). The prominence of this feature in Arp 2 20 is consistent with the presence of an OH 
megamaser in this source (ILonsdale et al.l 1 1 9941 ). 
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3.2. Unusual lines 



In addition to the lines discussed in the previous section, we also see a variety of 'rare' 
(which we arbitrarily define as appearing in ten or fewer objects) emission features in several 
objects, listed in Table HI We detect [C1II]A14.37 in six objects, four of w hich also show 
[NeV]A14.32, leading to a double-peaked profile (see also lArmus et al.ll2006l ). In the other 
two objects only one peak is seen, and it is possible that we have confused [C1II]A14.37 with 
[NeV]A14.32, however the velocity shift relative to the systemic (optical) redshift would have 
to be > 5000km s" 1 for this to be the case. We therefore regard [C1II] A14.37 as the more 



likely identification . This line is seen in so me lower luminosity starbursts (ISpoon et al.ll2000 



Verma et al.l 120031 1 . but is rare in AGN (ISturm et al.ll2002l ). A number of low-ionization 
iron lines are present, including [FeII]AAA17.94,24.52,25.99 and [FeIII]A22.93, all of which 
are seen (rarely) in lower luminosity starbursts and AGN. We tentatively detect [PIII]A17.89 
in two objects. This line is seen in small numbers of local IR-luminous sources, and in 
a variety of G alactic sources. A lso present is the [ArV]A13.10 line, which is occasionally 
seen in AGN ( ISturm et al.ll2002l ). One object, Arp220, shows a weak detection of what is 
plausibly [NeIII]A36.01, though Arp220 is the only object in our sample at a low enough 
redshift for this line to enter a well behaved part of the LH bandpass. We weakly detect 
the HI 7-6 line (i.e. the alpha transition of the Humphreys atomic Hydrogen series) in four 
objects. Finally, two objects show features that appear to be real, but proved difficult to 
identify reliably, hence we have not listed them in Table HI IRAS 23498+2423 shows a 
feature at a rest-frame wavelength of 10.581/im with a flux of 3.17xl0~ 22 W cm -2 . This 
feature is clearly resolved from the [SIV]A10.511 line. If this feature is real, then possibilities 
include |GoI I|A10.521, which i s seen in the 200 day - 500 day mid-IR spectra of core collapse 
supernovae faotak et alil2006h . or [Nill] A10.682. IRAS 19297-0406 shows an emission feature 
with a flux of 3.88xl0~ 21 W cm -2 that can plausibly be identified as the 31.77/xm water ice 
emission feature. 



4. Discussion 

4.1. Properties of the narrow-line region gas 

Mid-IR fine structure lines can be used to study three properties of the narrow-line 
region (NLR) gas§; excitation, electron temperature and electron density. Measuring electron 



3 Due to Spitzers limited spatial resolution, the 'NLR' of a ULIRG should be regarded as the ensemble of 
all the regions in a ULIRG where gas heated by star formation and/or an AGN emits fine-structure lines in 
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temperature though requires combining a mid-IR line with an optical or near-IR line of the 
same species. In this paper, we therefore concentrate only on electron density and excitation. 



4-1.1. Electron Density 

The electron density in a NLR can be estimated by using the ratios of two lines repre- 
senting transitions from the same orbital and suborbital (i.e. the same principal and angular 
momentum quantum numbers), but with different numbers of electrons in that (sub)orbital. 
For our sample there are two line ratios that satisfy this requirement; the 3 p2 — 3 pi and 
3 pi — 3 po transitions of [NeV] at 14.32/xm and 24.32/xm, and the 3 p2 — 3 pi and 3 p± — 3 po 
transitions of [SIII] at 18.71/xm and 33.48yum. These ratios can then be converted to an 
electron density by solving the relevant rate equation. Of the two, the [NeV] ratio is more 
useful. Its high ionization potential of 97eV means that it can only be produced (at a level 
observable in a ULIRG) by an AGN, making it more straightforward to interpret than the 
[SIII] ratio. Furthermore, for most e xtinction law s the extinctions of the two [NeV] lines are 



virtualy identical (Tabled see also iDraind 120031 ) . making the [NeV] ratio only marginally 



sensitive to the lack of an extinction correction. 

Sixteen of the ULIRGs in our sample show detections in both [NeV] lines. The line ratios 
for these objects span the range 0.57 < [NeV]A14.32/[NeV]A24.32 < 2.69, with a fairly even 
spread. The nine objects that are only detected in [NeV] A14.32 have lower l imits on the [NeV] 



line ratio consistent with this range. Using Figure 3 of lAlexander et all (119991 ). then these 
[NeV] line ratios are consistent with electron densities of < 10 4 cm -3 in all cases, well below 
the critical densities of the two lines. Our derived electron densities are also comparable t o 



those derived (using the same [NeV] ratio) for lower luminosity AGN (ISturm et al.l 120021 ). 
As a check, we examine the [SIII] line ratio; for those objects with detections in both lines 
we see a range of 0.14 < [SIII] A18.71 /[SIII] A33.48 < 1.10, consistent with electron densities 
of < 10 3,5 cm _3 . It is important to note however that we expect the [NeV] and [SIII] line 
ratios to give different electron densities, because the [NeV] lines are likely produced solely 
by the AGN, while the [SIII] lines (probably) arise from both starburst and AGN heated 
regions. 



the mid-IR 
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4-1.2. Excitation 



If electron densities are below the critical density, then the excitation level of the NLR 
can be estimated by considering flux ratios of adjacent ionization states of the same element, 
e.g. for an element X; fx^/fx*- F° r a fixed number of ionizing photons per Hydrogen 
atom (i.e. a given ionization parameter U), this ratio will be approximately proportional 
to the number of photons producing the observed X 1 flux relative to the number of Lyman 
continuum photons, though if U varies then fx^+i / fx* w iU a l so vary. Therefore, for a fixed 
U, a higher value of fx^ 1 / fx 1 indicates a harder radiation field, though the details depend 
on the element in question. 

Our spectra contain a limited number of lines, hence detailed excitation diagnostics 
are not possible. We can however get a qualitative idea of the range in excitations in 
our sample. The two most useful diagnostic ratios are [NeIII]A15.56/[NeII]A12.81 and 
[SIV]A10.51/[SIII]A18.71. Both Neon and Sulfur are abundant in Galactic sources, with 
both [Ne/H] and [S/H] lying approximately in the range 10 -4 — 10~ 6 . Furthermore, the pho- 
ton energies required to produce any of these four ions are < 50eV, meaning that they can 
all be produc ed in star-forming regio ns as well as A GN. All four lines ar e seen in Wolf-Rayet 
star spectra (ISmith fc Houckl 120011 ). HII regions (IPeeters et al.l 120021 ). planetary nebulae 
( Bernard Salas et al.ll200ll ). as well as the integrated spectra of local galaxies. We might 

SHIl in lower den - 
Ho fc Ketol[2007l ). 



also expect [Nelll] and [SIV] to increase in strength relative to [Nell] and 
sity and/or lower metallicity star forming regions (IThornley et al.l 12000 
It is however strange, given their similar ionization energies, that we detect [NeIII]A15.56 
in virtually all the sample, but only detect [SIV]A10.51 in about 45% of the sample. The 
[SIV]A10.51 line lies close to the 9.7/zm silicate feature, but we detect the H 2 ,S(3)A9.66 line 
in most of our sample, in many cases with a flux that is lower than the upper limit on the 
[SIV]A10.51 line. The most likely explanation for this is that the Neon and Sulfur emitting 
zones in ULIRGs lie wi thin regions extinct ed by silicate dust, while the H2 emitting region 
lies outside it (see also iHigdon et al.l 120061 ) . This is supporte d by the fact that we do not 
detect [SIV]A 10.51 in any o bject which has a silicate strength (ILevenson et al.ll2007l ) greater 



than 2.1 (see ISpoon et al.l (120071 ) for the silicate strengths of our sample) 



In Figure [TT1 we plot an 'excitation plane' of [NeIII]/[NeII] vs [SIV] /[SHI] for our sample. 
The ULIRGs span a broad range in excitation, but with a positive correlation between the 
Sulfur and Neon line ratios, suggesting that we are seeing emission from the same region 
in both s pecies. A simila r correlation has previously been noted f or Blue Compact Dwarf 
galaxies ( Wu et aD I2OO6I ) . and for nearby infrared-faint galaxies ( Dale et"aD I2OO6I ) . The 
ULIRGs are fairly evenly distributed across the correlation region, with no obvious over- or 
underpopulated regions. 
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Also plotted in Figure [TTJ are da ta for starbursts and AGN with IR lum inosities (mostly) 
between 1 x 10 10 L^ and 3 x 10 n Lp JVerma et alJl2003l : Isturm et al.ll2002f ). and the relations 
from figure 4 of iDale et al.l (120061 ) for star forming regions and IR-faint AGN, the bulk of 
which have IR lumino sities below 1 x 10 1 L ffi . The ULIRG s occu py the same region in the 
ionization plane as the lVerma et al.l (120031 ) and lSturm et al.l (120021 ) starburst and AGN sam- 
ples. Barring some particularly subtle bias from a lack of extinction correction, this suggests 
that the mid-IR emitting narrow line region gas in ULIRGs is excited in broadly the same 
way as the NLR in 10 10 < Lj r (L Q ) < 3 x 10 11 starbursts and AGN, and that any differences in 
density, metallicity or star formation history between ULIRGs and lower luminosity systems 
are not sufficient to manifest themselves in simple ionization plane diagrams. 



The same does not however appear to be true for sy s tems with L ir < 10 10 L 



The 

ULIRGs are consistent with the slopes of the IDale et al.l (120061 ) relations, but are offset 
below them by ~0.2 dex. Increased extinction in the form of a foreground screen could 
explain this, but would require an additional Ay ~ 60 of foreground extinction (and see also 
Table [2]). An (arguably) more likely explanation is that there is a decrease in the intrinsic 
[NeIII]/[NeII] ratio for a given [SIV]/[SIII] ratio in going from Lj r ~ 1O 1O L starbursts to 
L ir ~ 10 12 L Q systems. It is beyond the scope of this paper to investigate this effect in detail, 
so here we simply suggest a possible explanation. From Table EJ [SIV] has a significantly 
smaller ionization energy than [Nelll] , and the difference in ionization energies between [SIV] 
and [SIII] is smaller than the difference between [Nelll] and [Nell]. This means that, if we 
increase the gas density in the NLR or decrease the hardness of the ionizing radiation then 
both the [NeIII]/[NeII] and [SIV] /[SIII] ratios will decrease, but the [NeIII]/[NeII] ratio will 
decrease by a larger fraction than the [SIV] /[SIII] ratio. As we might expect the gas density 
in the NLR of ULIRGs to be significantly higher than in Lf r < 1O 1O L star forming regions, 
we might also expect a smaller [NeIII]/[NeII] ratio in ULIRGs for a given [SIV] /[SIII] ratio 
than in lower luminosity starbursts. A softer ionizing radiation field in ULIRGs compared 
to lower luminosity starbursts would produce a similar effect. 

We can also explore variations in excitation as a function of the total emission. In Figure 
[T2l we plot the [NeIII]/[NeII] ratio against infrared (rest-frame 1 — 1000/im) luminosity. 
There is no discernible correlation. Two ULIRGs, Mrk 463E and 3C 273, have elevated 
[NeIII]/[NeII] ratios (though with very different IR luminosities), while the rest of the sample 
is confined to [NeIII]/[NeII] ratios in the range 0.08 - 1.05, irrespective of their IR luminosity. 
This picture remains broadly the same if we substitute IR luminosity for radio luminosity 
(Figure [T3"j) . In this case, objects with 1.4GHz luminosities lower than ~ 10 24,5 W show no 
discernible correlation between radio luminosity and [NeIII]/[NeII] ratio. Those objects with 
higher radio luminosities appear to be confined to higher [NeIII]/[NeII] ratios, and with a 
narrower spread, though the number of objects with L1.4 > 10 24,5 W is too small to draw firm 
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conclusions. We conclude therefore that neither 1-1000/im luminosity or 1.4GHz luminosity 
are good proxies for the excitation of the NLR in ULIRGs. 



4.2. Starburst & AGN activity 



Mid-IR emission from galaxies can arise from five sources; (1) a non-thermal compo- 
nent, e.g. supernova remnants, (2) photospheres of evolved stars, (3) ionized gas, (4) dust 
grains, and (5) molecular gas. In most ULIRGs the latter three sources dominate the mid- 
IR emission. The ionized gas gives rise to the fine structure lines, while the smaller dust 
grains produce the mid-IR continuum longward of about lOyum. La rger dust grains give rise 



to features such as the absorption features at 9.7/im and 18/xm (ILebofsky fc Riekd Il979 



Roche fc Aitkenl Il985l ; IChiar fc Tielend 120061 ) . Large molecules give rise to several emission 
and absorption features; the most prominent emission features are seen from 3.3/xm to about 
19/zm and arise from bending and stret ching modes of Polycyclic Aromatic Hydrocarbons 
( jPuget fc Legerl Il989l ; iHony et al.l l200ll ) . Other molecular features include several Hydro- 
gen emission lines, as well as weaker feature s such as CO, HCN and C2H2 gas absorption 
( ISpoon et al.ll2006l ; iBernard-Salas et al.ll2006l ; lLahuis et al.l 120071 ). 



The use of these mid-IR emission and absorption features as diagnostics of the power 
source in ULIRGs is however not straightforward. The fundamental problem is that dust and 
gas simply require a source of ionizing radiation to emit in the mid-IR, and do not particularly 
care what the source of that ionizing radiation is. In principle, a population of hot stars 
or an accretion disk around a supermassive black h ole both serve equally well. Exhaustiv e 
reviews of these difficulties are given elsewhere (e.g. lLonsdale et al.ll2006l ; iDale et al.l 120061 ) . 
so here we highlight some examples. Fine structure line fluxes are affected by several factors, 
particularly in starbursts, where the age of the burst, the geometry of the region and the 
upper and l ower bounds on the I MF can all have a drastic effect on observed line ratios 
(e.g. fig 6 of iThornley et al.ll2000l ). The factors determining PAH strengths are still poorly 
understood; one factor that is known to have a significant effect is metallicity, with strong 



supp r ession of PAHs seen in low metallicity environments (IRoche et al.l 1 199 It iThuan et al. 



1999; 


Dwek 


2005; 


Wu et al. 


2006) 



stretching modes cannot be excited by an AGN, although observationally PAHs seem to be 
almost exclusively a ssociated with either star forming regions or ambient interstellar light 



( iPeeters et al.l 120041 ) . Significant silicate absorption on the other hand merely requires a 
mass of warm dust obscured by a significant column of cooler dust, and therefore on its own 
says nothing about what is heating the warm dust. 



With these subleties in mind, we examine the nature of the starburst and AGN activity 
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in our sample. We focus here on a limited number of diagnostics, deferring detailed analysis 
to upcoming papers. 



4-2.1. Neon lines and star formation 



Recently, iHo fc Ketol (120071 ) have suggested that the total luminosity of the [NeII]A12.81 
and [NeIII]A15.56 lines are a diagnostic of star formation rates in galaxies across a wide range 
in IR luminosity. This method potentially offers a number of advantages over other mid-IR 
star formation diagnostics; the Neon lines in question are strong and easily observed, and 
the dependence on an accurate extinction correction is small. With ULIRGs however, there 
is the obvious caveat that some or all of the [NeII]A12.81 and [NeIII]A15.56 may originate in 
regions excited by an AGN rather than a starburst, so in this section we explore the use of 
the [NeII]A12.81 and [NeIII]A15.56 lines as star formation rate diagnostics in our sample. 

In Figure we plot the [NeIII] + [NeII] luminosities of our samp l e aga inst total IR 
luminosity, and overplot the relation given in equation 2 of IHo &: Ketol (120071 ). There is a 
large scatter in the [NeIII] + [NeII] luminositie s as a functi o n of I R luminosity for our sample, 
comparable to the scatter seen in figure 1 of IHo fc Ketol (120071 ). but we see a clear upward 
trend in [NeIII] + [NeII] luminosity with increasing IR luminosity. A power law fit yields 
Ltv oc L^ 75 , and a horizontal line fit (i.e. no dependence of [NeIII] + [NeII] luminosity on IR 
luminosity) is ruled out at > 5cr significance. The slope of the relation for the ULIRGs is 
consistent with the slope of the Ho & Keto relation plotted in Figure H~4"l but the ULIRGs are 
systematically lower, offset by ~ 0.4 dex, on average. There are several possible explanations 
for this offset. First is that an AGN is contributing to the total IR luminosity but not to 
the Neon line fluxes, however this explanation seems unlikely, as we see systems with known 
IR-luminous AGN scattered on either side of the Ho & Keto line; Mrk 463E and NGC 6240 
both l ie above it, whereas Mrk 231 and IRAS 03158+4227 (which contains a Compton-thick 
AGN, iRisaliti et al.ll2000l ) both lie well below it. Second is a different electron density in the 
NLR of ULIRGs; either significantly lower than that seen in lower luminosity starbursts, or 
much higher and approaching the critical density. A high electron density however is ruled 
out (see §4.1.11) . and there is no evidence for electron densities substantially lower than in 
lower luminosity systems. A significantly lower metallicity would also serve to lower Neon 
luminosities for a given IR luminosity, but again there is no evidence (or indeed plausible 
motivation) for such low metallicities in local ULIRGs. Finally, this offset could be caused by 
higher extinction. This explanation seems the most likely. Though the Ho & Keto relation 
was derived without correcting for extinction, we expect ULIRGs to suffer heavier extinction 
in their nuclear regions relative to lower luminosity starbursts. The offset seen for ULIRGs 
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corresponds to Ay = 40 ± 20 over and above the extinction seen in lower luminosity systems. 



4-2.2. PAHs and star formation 



The correlation between [NeIII]A15.56 + [NeII]A12.81 luminosity and IR luminosity is 
consistent with the idea that there is a correlation between [NeIII]A15.56 + [NeII]A12.81 
luminosity and star formation rate in ULIRGs. This on its own however does not support 
such a hypothesis, as we expect the luminosity of any mid-IR fine-structure line to correlate 
approximately with the mid-IR luminosity of the ionizing souce, irrespective of what that 
source may be. Therefore, to test this hypothesis further, we examine the relationship 
between [NeIII]A15.56 + [NeII]A12.81 luminosity and the luminosity of PAH features. 



The origin of PAH features and their use as star forma tion indicators are still controver 



sial 



e.g. 



Peeters et al.ll2004l ; iForster Schreiber et al.l I2004T ) . Observationally however, PAH 



features are usually prominent in starburst galaxies, with previous authors noti ng a good cor- 



relati on between the strength of PAH features and IR luminosity in starb ursts (IBrandl et al. 



20061 ). but appear to be weak or absent in AGN (jWeedman et al.ll2005l ). The cause of this 
dichotomy is thought to be a combination of two factors. First, the UV radiation field from 
an accretion disk around a black hole is harder than the UV radiation field from a starburst, 
and is therefore more adept at destroying the C-H and C-C bonds in PAHs. Second, a lu- 
minous AGN produces a prominent mid-IR continuum which ca n drown out PAH features, 
even if there is a vigorous starburst present (ILaurent et al.l 120001 ) . 



Even though there is good evidence that PAHs originate in star forming regions, the 
behaviour of individual PAH features as a function of star formation rate remains uncertain, 
with strong var iations observed in individual PAH feature strengths between different star- 
burst galaxies (jSmith et al.l 120071 ) . Therefore, we elect to use the combined luminosities of 
two PAH features in order to reduce the likely scatter introduced by variances in individual 
PAH strengths. The choice of the two PAH features is a straightforward one; the most lumi- 
nous PAH features are those at 6.2/im, 7.7/im and 11.2/im, but the 7.7/iin feature is difficult 
to measure as it lies in a crowded part of the spectrum. Therefore, even though the 11.2/im 
feature can be significantly affected by silicate absorption at 9.7/im, we use the luminosities 
of the PAH 6.2/im and 11.2/im featu res. These PAH lum inosities are measured from the low 
resolution IRS spectra presented in iDesai et al.l (120071 ). For the 11.2/im PAH feature the 
luminosity was measured by integrating over 10.8/im - 11.8/xm in the continuum-subtracted 
spectra, with the equivalent width evaluated at 11.25/xm. The relevant parameters for the 
6.2/im PAH feature were 5.90/im - 6.55/im and 6.22/im, respectively. 



-16- 



In Figure [15] we plot [NeIII] + [NeII] luminosity against the combined luminosities of the 
PAH 6.2/im and 11.2/xm features^. The Neon line luminosities and PAH luminosities clearly 
track each other, lending support to the idea that, on average, lower ionization Neon lines 
originate in regions heated mostly by star formation. A power law fit yields: 



(0.17 



+0.46\ 
0.12 J 



x L 



1.02±0.05 



(3) 



where and Lp are the [NeIII] + [NeII] and 6.2/zm+11.2 / um luminosities (in any units, as 
long as they are both the same), respectively. Equation [3] is solely intended to indicate the 
trend, and evidently is not a good physical model for the data. It is notable however that 
the scaling between [NeIII] + [NeII] luminosity and PAH 6.2/im+11.2/zm luminosity predicted 
from this fit is, to within the error on the exponent of Lp, linear. Interestingly, we see a 
comparable slope (though different normalization) if we instead plot [SIII] A18.713 luminosity 
against PAH 6.2/im+11.2yum luminosity. This is consistent with the idea that the sizes of 
the lower ionization Neon and Sulfur emiting regions, and PAH emitting regions in ULIRGs 
scale linearly with each other with increasing luminosity, though a luminosity dependence 
on metallicity and/or electron density could also play a role. 

It is also straightforward to derive a relati on between PAH 6 .2/xm+11.2/xm luminosity 
and star formation rate, using equation 12 from lHo fc Ketol (120071 ) combined with Equation 
El In doing so however, we u se the updated convers ion between Lyman continuu m flux and 
star fo rmation rate given by iHirashita et all (120031 ). rather than the relation in iKennicutt 



( 1998 ). By combining equation 12 from Ho Sz Ketol ( 2007) with the relat i ons be tween Lyman 
continuum flux and star formation rate given in §2.1 of IHirashita et al.l (120031 ). we obtain: 



SFR[M Q yr- 1 } = 2.69 x 10~ 41 



L at [ergs s~ 



fion(fNe+ + l-67/jVe++) 



(4) 



where SFR is the star formation rate, L^ eon is the combined luminosity of the two Neon lines, 
f ion is the fraction of photons that actually contribute to ionizing the gas, and f Ne + and fwe++ 
are the fractional abundances of [Nell] and [Nel H], respectively. For U LIRGs however there 
is an additional complication. The conversion in IHirashita et al.l (120031 ) assumes continuous 
star formation over > 10 8 years. Star formation in ULIRGs however is likely to occur in some 
form of 'burst' over timescales of 10 7 — 10 s years. This means that the star formation rates 
in ULIRGs will be underestimated by a significant amount if the conversion in Equation H] is 



4 A similar plot but using only the 6.2^m luminosity or the 11.2/xm luminosity yield plots with significantly 
greater scatter 
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used unmodified. Following previous work (jKennicuttlll998l ; iThornley et al.ll2000l ). we adopt 
an approximate upward scaling of 50% to correct for this difference. This scaling is likely to 
be sufficient for considering trends, but should be adjusted for specific objects where the age 
and nature of the burst are known. Assuming this scaling, plus linear scaling bet ween Neon 
and PAH luminosity, fi on = 0.6, f^ e + = 0.75 and fNe++ = 0.1 (IHo fc Ketol 120071 ). we arrive 
at: 



SFR[M Q yr~ 



1.18 x l0~ AL L P [ergs s~ 



(5) 



where Lp is in units of ergs s _1 . This assumes Solar metallicity and a Salpeter IMF spanning 
O.1M to 100M o , and is only applicable to objects where the PAH and Neon luminosities 
arise from recent star formation. 

The errors on the star formation rates derived using Equation [5] are difficult to quan- 
tify. The convers ion from Lya continuum flux to star formation rate has an error of ~30% 
( IKennicuttl Il998l ) . and there is a significant extra error introduced by the scatter in Figure 
IT51 We therefore estimate that the errors on star formation rates derived using Equation [5] 
are of order 30% for population studies, rising to at least 50% for individual objects. Equa- 
tion [5] does however offer the advantage that it does not require an accurate measure of the 
rest-frame 1-1000/xm luminosity, and hence is particularly suited to population studies at 
high redshifts. Equation [5] also gives plausible star formation rates for our sample. We defer 
a complete ana lysis to an upcomi ng paper, but note that for Arp 220 and using the PAH 
luminosities in iDesai et all (120071 ). Equation [5] gives a star formation rate of 57M Q yr _1 , 
which is consistent w ith the total star f ormation rat e inferred from rad io observations, which 
give 5O-1OOM yr" 1 rtSmith et allll998l . but see also lParra e^al1l2006h . 



4-2.3. Starburst vs. AGN diagnostics 

To distinguish between starburst and AGN power behind the 1-1000/im emission in 
ULIRGs, we need diagnostics that are sensitive to the observable differences between young 
stars, and an accretion disk around a supermassive black hole. For our purposes, there are 
two such differences. First, the ionizing radiation from an AGN is harder and (potentially) 
more intense than that from a starburst. Second, AGN occupy a smaller volume, ~ O.lpc, 
as opposed to a few tens to a few hundreds of pc for a starburst. In principle therefore, 
we might expect the averaged mid-IR spectra of AGN to exhibit three differences compared 
to those of starbursts; higher ionization fine structure lines, increased quantities of 'hot' 
(> 100K) dust, and (potentially) increased obscuration. We explore each of these differences 
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in this section. 



We start with the simplest possible diagnostic; the detection (or otherwise) of individual 
lines. For examining starburst and AGN activity, the [NeV]A14.32 and [OIV]A25.89 lines 
are the most useful single line di agnostics. Both l ines can be strong in planetary nebulae 
and young supernova remnants (jOliva et al.lll999l ). and [OIV] is sometimes se en in W-R 



star sp ectra. Both lines are however weak in spectra of star forming regions (e.g. iLutz et al. 
( 11998I )). while being strong in spectra of AGN. As [NeV]A14.32 has E ion = 97.1 eV, compared 



to 54.9eV for [OIV]A25.89, and becau se models i ndicate that [NeV]A14.32 is unlikely to be 
strong in galaxies without an AGN fjVoitl Il992l ). we focus most of our attention on the 
[NeV]A14.32 line. 

The presence of [NeV]A14.32 cannot provide meaningful constraints on the bolometric 
luminosity of an AGN, simply becase a comparable [NeV]A14.32 flux can be produced by a 
faint AGN with low obscuration, or a luminous AGN with high obscuration. From Table [3] 
however we can see that [NeV]A14.32 is detected in 22 out of 53 objects, providing direct 
spectroscopic evidence for an AGN in these 22 objects. Our sample is flux- rather than 
volume-limited, hence we cannot draw firm conclusions about local ULIRGs as a whole, 
however this suggests that an AGN provides a non- negligible fraction of the mid-IR flux 
in ~ 42% of local ULIRGs. This is significantly higher than the fraction of ULIRGs that 
show spectroscopic signatures of AG N activity in their optical o r near-IR spectra, 20% — 
25% (IVeilleux. Sanders &: Kimlll999i ). though strictly speaking the IVeilleux. Sanders fc Kim 
19991 sample and our sample are not directly comparable as their selections are different. 
Interestingly, we see no convincing trend in the detection of [NeV]A14.32 as a function of 
IR luminosity. From Table [3] the ULIRGs with a [NeV]A14.32 detection are fairly evenly 
spread between IR luminosities of 1O 1L8O L and 1O 1288 L ; 10/22 ULIRGs with a [NeV]A14.32 
detection have L ir > 1O 12,38 L , while 20/53 of the whole sample have L ir > 1O 12 38 L . 
This contrasts with optical/near-IR spectroscopic surveys, which report a rising fraction of 
ULIR Gs with AGN signatures, reaching 35% -50% at > 1O 12 38 L (IVeilleux. Sanders fc Kim 
A plausible reason for this discrepancy is high levels of extinction toward the nuclei 
of ULIRGs. From §4.2.1[ the increase in extinction toward the nuclei of ULIRGs compared to 
lower luminosity systems is of the order Ay ~ 40, making it harder to see AGN signatures in 
the optical or near-IR. Mid-infrared spectroscopy on the other hand appears to be a sensitive 
probe of the presence or otherwise of AGN in IR-luminous systems. We detect [OIV]A25.89 
in 21/22 objects that show [NeV]A14.32, and detect [OIV]A25.89 in only two objects that 
do not show [NeV]A14.32, suggesting that [OIV]A25.89 is a good, but not perfect proxy for 



5 the 'boundary' luminosity quoted by IVeilleux. Sanders fc Kind 11999 is 1O 12 3O L0, however this is for 



H = 75 km s _1 Mpc -1 and q = 0.0 
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[NeV]A14.32. 



We next consider diagnostics based on fine-structure line ratios. In Figures [H] and [T7] 
we plot [NeV]/[NeII] and [OIV]/[NeII] against IR luminosity, together with the predicte d 
AGN contribution to the total IR luminosity based on these line ratios (ISturm et al.ll2002l ). 
Both diagrams predict a broad spread in power source, ranging from 100% starburst to 100% 
AGN. In both diagrams however only 10-12 ULIRGs are above the 40% AGN line, and more 
than half the sample lie below the 20% lines. This is consistent with star formation being 
the dominant contributor to the IR emission in most ULIRGs, with only 20% of ULIRGs 
hosting an AGN with a comparable or greater IR luminosity than the starburst, in agree- 



ment with studies at ot her wavelengths (jVeilleux. Sanders fc Kiml Il999i ; iFarrah et all 12003 



Franceschini et all 120031) . There are however many caveats in using such simple diagnostics 
(ISturm et al.ll2002l ). and Figures [TBI and [T71 should not be considered reliable diagnostics 
of the power source in individual ULIRGs, but rather as crude indicators of trends. It is 
however notable that we seem to have results consistent with other diagnostics for (some) 
individual objects. For example, both diagrams predi ct a small contr i butio n from an AGN 
i n Arp 220, consis tent with results from the X-ray (iClements et al.l l2002l ) and from ISO 
Jsturm et alJll996h . and 3C 273, Mrk 1014, Mrk 463E and IR AS 05189-2524 are all pre- 



dicted to contain a luminous AGN, in line with previous w ork (jShakura fc Sunyaevl Il976 
Mazzarella et alill99ll : iBoller et al1l2002l : IFarrah et allboosh . 



Given the likely origin of the [NeV]A14.32 line in AGN, we can use the [NeV]/[NeII] 
ratio to test other mid-IR AGN diagnostics. One such diagnostic is the IRAS 25/im/60/zm 
flux density ratio; it has been suggested that 'warm' objects (those wi th / 2 n/ 'fen > 0.2) are 
more likely to con tain an AGN than 'cool' (/25//60 < 0.2) objects (jde Grijp et al.l Il985l ; 
Sanders et allll988l ). In Figure [T8l we plot [NeV]/[NeII] against fzzl fm- Just over half of 



the 'warm' objects have [NeV]A14.32 detections, compared to about one third of the 'cool' 
objects. Furthermore, for those objects with [NeV]A14.32 detections, the 'warm' objects have 
systematically higher [NeV]/[NeII] ratios than the cool objects. While our small sample sizes 
(particularly for the warm objects) render any conclusions tentative, we infer from this that 
the IRAS 25/im/60/im flux density ratio is a reasonably good, though not perfect diagnostic 
for the presence of an IR- luminous AGN in a ULIRG. 

We move on to consider diagrams with a diagnostic of the starburst and AGN lumi- 
nosities on both the x and y axes. A natural diagnostic to combine with a fine-structure 
line ratio is a measure of the strength of the PAH features. To measu re PAH strength w e 
use the equivalent width (EW) of the 6.2/zm PAH feature, taken from iDesai et al.l (120071 ). 
as this feature is strong, and lies in a relatively uncluttered part of the mid-IR spectra of 
ULIRGs, though its use does depend on an accurate correction for water ice and/or aliphatic 
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hydrocarbon absorption. A mid-IR fine-structure line ratio vs. PAH EW diagno stic diagram 
has been used by severa l authors for many classes of IR-luminous galaxy (e.g. 
19981 : Isturm et al.ll2002h . 



Genzel et al 



In Figures [Hand [20] we plot [NeV]/[NeII] and [OIV]/[NeII] against the equivalent width 
of t he 6.2/im PAH feat ure, together with the linear mixing ratios from ISturm et al.l (120021 ) 
and lArmus et al.l (120071 ). Here we again see a broad spread in fractional AGN and starburst 
luminosities, from 'pure' AGN to 'pure' starbursts. Those sources classified by Figures [Inland 
[2U] as AGN dominated tend to have PAH|^ < 0.05, [NeV]/[NeII]> 0.2 and [OIV]/[NeII]> 
0.5, whereas the 'pure' starbursts tend to have PAHj^ > 0.2, [NeV]/[NeII]< 0.09 and 
[OIV]/[NeII]< 0.2. The IR emission from the majority of the sample are still predicted to 
be powered mainly by star formation, though this trend is not so clear as from Figure [HJ 
We see broad (to within a factor of three) agreement between the two axes, and also that no 
object lies in the 'forbidden' top right hand corner of either diagram. In both diagrams, the 
six objects which show a detection of the [C1II]A14.37 line ( §3.2p are mostly located toward 
the right hand side, consistent with the postulated origin of this line in star forming regions. 

It is however notable that several sources lie in a 'forbidden' region in each diagram, 
where the PAH equivalent widths suggest a low starburst contribution, but the fine-structure 
line ratios suggest a low AGN contribution. There are, broadly, four possible reasons for 
this discrepancy^; (1) underestimated 6.2/xm PAH EWs due to nearby water ice and/or 
aliphat ic hydrocarbon absorption, (2) destruction of PAHs in lum inous starburst environ- 
ments (jRigopoulou et al 1999 ; Farrah et al. 2003 ; Beirao et al. 20061 ). (3) suppression of PAH 
EWs in compact starburst environments, and (4) increased general obscuration levels. We do 
not consider the first two possibilities likely though. We have in all cases used ice absorption- 
corrected values for the 6.2/im PAH EWs (the correction is o f order 10% for the 6.2/zm feature 
and less than 5% for the 11.2/im feature. ISpoon et al.ll2007l ) and it is unlikely that any resid- 
ual ice absorption could shift sources by such a significant distance on this plot. If the PAHs 
were being destroyed in extremely luminous starbursts then we would expect the objects in 
the 'forbidden' region to all have high IR luminosities, but the objects in this region (e.g. 
IRAS 08572+3915, IRAS 00397-1312, IRAS 15462-0450) span a wide range in IR luminos- 
ity. A combination of increased overall obscuration and increased starburst compactness 
therefore seem the most likely culprits. The mixing ratios in both diagrams are based on 
starbursts and AGN with IR luminosities of 10 10 — 10 11 L Q . Their use therefore assumes a 
linear scaling between starburst and AGN luminosity, and overall obscuration and starburst 



6 A fifth possibility is variation in the filling factor of the coronal-line region clouds. A low filling factor 
would imply a small ionization parameter U , and would therefore give an unusually low [NeV]/[NeII] ratio, 
for example. We do not however have the data to address this possibility 
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geometry, to give rise to comparable observed [NeV]/[NeII] and [OIV]/[NeII] ratios and PAH 
6.2/im EWs. These assumptions are unlikely to be valid. For exam ple, previous authors have 
noted A\'s of ten or more, even in the mid-IR (IGenzel et allll998l ). and we derived Ay ~ 40 
over and above lower luminosity stabursts in §4.2.11 Increased overall obscuration will lead 
to greater suppression of the [NeV] and [OIV] lines compared to the [Nell] line, as these 
regions must lie closer to the starburst and/or AGN. This increased obscuration will also 
arise in lower than expected PAH luminosities, but is unlikely to affect their EWs. However, 
we expect more compact starbursts to have a stronger mid-IR continuum, which will lead 
to suppression of PAH EWs. Overall therefore, increased overall obscuration combined with 
more spatially compact starbursts in ULIRGs compared to lower luminosity systems can 
plausibly explain the outliers in Figures [19] and I 



An important caveat however is that diagrams like Figures dH] and [2D] are not suited 
to studying the physics of individual objects, and can often get the answers wrong. For 
example, if the power source is so obscured that it does not emit significantly at wavelengths 
shortward of ~ 25/im, t hen mid-IR fine stru cture line ratio vs. PAH EW diagnostic plots can 
give misleading results (jPeeters et al.ll2004l ). Some 'outliers' are therefore understandable. 
One example is IRAS 15206+3342 (#28). This object lies squarely in the 'starburst' part 
of both plots, and from this one m ight conclude that it does not contain an energetically 
significant AGN. Its UV spectrum (IFarrah et al.l 120051 ) however shows clear evidence for a 
broad absorption line QSO. A second example is NGC 6240, which also lies in the 'starburst' 
part of both plots, despite the extensive evidence for an IR-luminous AGN in this source 
(lArmus et al.ll2006l ). The key point is that diagnostics such as those in Figures [T91 and |2"U1 
probe the physics of specific environments within ULIRGs. To obtain a complete picture of 
an individual object therefore requires multiple diagnostics to probe multiple environments. 



4-2.4- Silicate absorption 

We explore the issue of mid-IR obscuration in ULIRGs further by combining the diag- 
nostics discussed previously with a measure of the strength of the 9.7/im silicate absorption 
feature, S s n: 

= ln ( F obs (9.7nm) \ 
\F cont (9.7/j,m) J 

where F f> s is the observed flux at rest-frame 9.7/zm, and F cont is the underlying continuum 
flux at rest-frame 9.7/im deduced from a spline fit to the continuum flux at rest frame 5.0-7.0, 
14.0-14.5, and 25.0-40.0/im. A complete description of the method used to measure S si i can 
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be found in lSpoon et al.ll2006l and lLevenson et al.ll2007l . The silicate strengths for our sa mple 
are measured from the low resolution spectra, and are presented in ISpoon et al.l (120071 ). 



In Figure I2T1 we plot S^u against the [OIVl/[NeII1 line ra tio, along with the 'reference' 
starbursts and AGN farandl et al.lbood : IWeedman et alJ2005h . The reference starbursts and 
AGN separate well on the y axis, with most of the AGN above [OIV]/[NeII]=0.2 and most of 
the starbursts below. Nearly all the starbursts and AGN are confined to —0.5 < S s u < 2.0. 
The ULIRGs are found over the whole range of [OIV]/[NeII] ratios where the reference 
samples are seen, but are offset on the x axis, with greater silicate strengths for a given 
[OIV]/[NeII] ratio. The simplest interpretation of this is that the starbursts in ULIRGs are 
similar in nature to those in lower luminosity systems, but with significantly greater (Ay's 
of a few tens) total obscuration, consistent with the conclusion in §4.2.11 

It would be natural, from this, to postulate that silicate strength is a measure of the 
total obscuration of the starburst. If so, then we might expect a correlation between star 
formation rate and silicate strength; the reasoning being that a more luminous starburst 
will have a deeper 9.7/zm absorption feature as the starburst has an increased total dust 
column. It is straightforward to test this hypothesis; in Figure [22] we duplicate Figure [Jj)] 
but this time color-coded each object according to its silicate strength. Clearly, there is no 
correlation. Objects with the strongest silicate absorption seem to be confined to the lower 
left of the plot, though with only three objects it is impossible to draw firm conclusions. 
The other objects seem to be scattered randomly; strongly absorbed and weakly absorbed 
sources are found from the top right to the bottom left. Furthermore, if S s u was simply 
a measure of the obscuration of the starburst, then we would expect the sources with the 
greatest negative offset from the Ho & Keto line in Figure [15] to have the largest values of 
S si i, but this is not the case. For example, three of the sources with the greatest negative 
offset in Figure [151 IRAS 06035-7102, IRAS 20551-4250, and 3C 273, range from having a 
deep silicate absorption feature to a silicate emission feature. We infer from this that silicate 
strength is not simply a measure of the total ob s curat ion of the starbursts in ULIRGs. A 
similar conclusion was reached by iHigdon et all (120061 ) , based on the molecular Hydrogen 
line strengths. 

We move on to consider a different scenario, in which the 9.7/im silicate absorption in 
ULIRGs is affected by the AGN as well as, or instead of, the starburst. To examine this 
possibility we duplicate the plot in Figure [20] in Figure [23] again color-coding each point 
by its silicate strength. From this we see that both strongly and weakly absorbed sources 
are found across the span of the plot, and there is no appreciable trend in average silicate 
strength in any direction. There is however one interesting trend. Those sources towards 
the left side of the diagram tend to have either very deep absorption, or very shallow to 
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no absorption, whereas those objects on the right hand side of the plot seem to be almost 
universally 'averagely' absorbed. Of the 13 objects on the right hand side, all but two (IRAS 
15206+3342 & IRAS 23128-5919) have 0.8 < S sll < 2.4. Of the 13 objects on the left hand 
side, only one (IRAS 12514+1027) does not satisfy S s u > 2.4 or S s u < 0.8. Overall, there is 
a clear shift in the distribution of values of S s u as we move from the left hand side to the 
right hand side of the plot, from a single, broad peak centered at around S s u — 1.5 to two 
narrower peaks, one at S s u ~ 2.8 a nd on e at S a u ~ 0.2. This shift is conceptually similar to 



the 'fork' diagram in ISpoon et al.l (120061 ). 



These distributions in silicate strength can be interpreted in one of two ways. The 
first is that very strong or very weak silicate absorption indicates the presence of an AGN, 
with moderate silicate absorption indicating the presence of a starburst. For an AGN we 
expect a strong dependence of observed properties on viewing angle, as nearly all lines of 
evidence suggest tha t the dust in AGN is a rranged in a pl anar structure, w hether that 
structu re is a torus (ISchartmann et al.l 120051 ) . a flared disk (IFritz et al.l 120061 ). or discrete 
clouds (jElitzur fc Shlosmanll2006l ). For a starburst however we expect a weaker (though still 
possibly significant) dependence on viewing angle; millime ter interferometry h as shown that 
the starburst re gions in UL IRGs are dense and compact (ITacconi et al.l 120061 ) . while recent 



theoretical (e.g. ISilk! (120051 )) and indirect observati onal eyidenc e suggests that starbursts in 



ULIRGs are unlikely to have a disklike structure (jFisherl 120061 ) . We therefore expect that 
the distribution of values of S s u in obscured AGN will be bimodal, with a high apparent 
obscuration when the AGN is viewed edge on, but low obscuration when the AGN is viewed 
face on, but that the dependence of S s u on the viewing angle of the starburst will be weaker 



This is supported by the fact that optical QSOs usually show silicates in emission (IHao et al. 



20061 ). and that one of the two objects on the right hand side of the lower panel of Figure 
is IRAS 15206+3342 which, as previously mentioned, contains a BAL QSO. 



The second way is that the sources with S s u > 2.4 instead contain an extremely obscured 
starburst, and that only those sources with either a [NeV] detection or a silicate feature in 
emission contain an IR- luminous AGN. This scenario is plausible as a very compact, highly 
obscured starburst could also result in suppression of PAH features, and would also explain 
the absence of [NeV] detections (if the starburst surrounds the AGN) in the few objects on 
the left hand side of Figure [2"31 with S s u > 2.4. With the data available to us it is difficult to 
ch oose between these two possibilities. Both scenarios are consistent with the 'fork' diagram 



m 



Spoon et al.l (120061 ) for example. With some reserve therefore, we propose that the first 



scenario is more likely, and that moderate silicate absorption signifies the presence of a 
dominant starburst, but that very deep or very shallow silicate absorption in ULIRGs likely 
signifies the presence of a bolometrically significant AGN. 
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5. Conclusions 

We have presented an atlas of fine structure lines and other emission features measured 
from high resolution mid-infrared spectra of 53 ULIRGs at 0.018< z <0.319, taken using 
the Infrared Spectrograph onboard the Spitzer space telescope. We have employed a variety 
of diagnostics using these emission lines as well as those based on PAH features and the 
strength of the 9.7/zm silicate absorption feature to investigate the power source behind the 
infrared emission. Our conclusions are: 

1) All of the spectra show various of fine structure emission lines of Neon, Oxygen, and Sul- 
fur, as well as one or more molecular Hydrogen lines. We see the [NeII]A12.81, [NeIII]A15.56 
and [SIII] A18.71 in most of the sample, and [SIV]A10.51 in just under half the sample. 
The higher ionization lines [NeV]A14.32 and [OIV]A25.89 are detected in just under half 
the sample, while [NeV]A24.32 is detected in about one third of the sample. Rarer lines 
include [C1II]A14.37 (six objects), HI 7-6A12.37 (four objects), [PIII]A17.89 (two objects), 
and [ArV]A13.10 (two objects). Some objects also show low-ionization iron lines, including 
[FeII]AAA17.94,24.52,25.99 and [FeIII]A22.93. The detection of three further lines is depen- 
dent on redshift, but we see [SIII] A33.48 in ~80% of the objects where this line lies within 
the bandpass, [ArII]A8.99 in just over 50%, and [SiII]A34.82 in just over 50%. 

2) The presence of the [NeV]A14.32/zm line in 22/53 objects is direct spectroscopic evidence 
for the presence of an AGN that provides a significant, though not necessarily dominant 
fraction of the mid-IR flux in ~ 42% of ULIRGs. Based on this, we find that the IRAS 
25 / um/60/im flux density ratio is a reasonable, though not perfect, diagnostic for the presence 
of an IR-luminous AGN in ULIRGs. In most, but not all objects where we see [NeV]A14.32/iin 
we also see [OIV]A25.89, suggesting that [OIV]A25.89 is a good proxy for the [NeV]A14.32/im. 
In contrast, we see [SIV]A10.51 in a surprisingly low fraction of the sample, given that its 
ionization energy is comparable to that of [NeIII]A15.56, and that the H2,S(3)A9.66 line 
is seen in most of our sample. The most likely reason for this is that the the Neon and 
Sulfur emitting zones in our samples lie within regions that are more strongly extincted 
by silicate dust than the H 2 emitting regions; the increased extinction due to the nearby 
9.7/xm absorption feature would then weaken the apparent flux of [SIV]A10.51 relative to 
[NeIII]A15.56. 

3) We use the [NeIII]/[NeII] vs [SIV]/[SIII] plane to show that the excitation levels in the 
mid-IR emitting regions span more than two orders of magnitude in both the Neon and 
Sulfur line ratios. The range in both line ratios is comparable to that seen in starbursts and 
AGN with IR luminosities in the range 10 10 <L ir (L ) < 10 115 , but we see a systematically 
lower [NeIII]/[NeII] ratio for a given [SIV]/[SIII] in our sample compared to systems with IR 
luminosities of < 10 10 L Q , possibly due to the increased density of gas in the NLR of ULIRGs. 
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We use the [NeV]14.32/im/[NeV]24.32/im and the [SIII]18.71/xm/[SIII]33.48/xm line ratios to 
show that the electron densities in the mid-IR emitting regions of ULIRGs are < 10 4 cm -3 
in all cases, well below the critical densities. 

4) We show that the combined luminosity of the [NeIII]A15.56 and [NeII]A12.81 lines cor- 
relates with both total IR luminosity (Figure [141) . and the combined luminosity of the 
PAH 6.2jnm and 11.2jLtm features (Figure fT5l) . Combining this result with previous work 



( ]Ho fc Ketd 120071 ). we derive a calibration between star formation rate and PAH 6.2/im + 



11.2/zm luminosity for ULIRGs: 



SFR[M G yr^] = 1.18 x 10~ 41 L P [ergs s' 1 ] (7) 

where M is the star formation rate in solar masses per year, and Lp is the PAH 6.2/mi + 
11.2/xm luminosity in ergs s _1 

5) We employ a variety of spectral diagnostics to show that, despite the presence of a 
luminous AGN in ~42% of ULIRGs, the most likely dominant contributor to the total IR 
emission in most ULIRGs is star formation, with an AGN providing a higher contribution 
than a starburst in only ~ 20% of ULIRGs. The fine structure line ratios, luminosities 
and PAH EWs of our sample are consistent with the starbursts and AGN in ULIRGs being 
more extincted (Ay ~ 40), and for the starbursts more compact, versions of those in lower 
luminosity systems. 

6) We show that the strength of the 9.7/zm silicate feature is unlikely to be a simple indicator 
of the total obscuration of the starburst. We combine measurements of PAH equivalent 
widths, 9.7/im silicate feature strengths, and fine structure line ratios to show that ULIRGs 
with silicate strengths of S s u < 0.8 likely contain an energetically significant AGN, whereas 
the IR emission from ULIRGs with 0.8 < S s u < 2.4 is likely dominated by star formation. We 
postulate that ULIRGs with S s u > 2.4 contain an deeply buried AGN, though a comparably 
obscured starburst is also possible. 
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Table 1. Observations summary 





Galaxy 


RA (J2000) 
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23 
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11 
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02 
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23 
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.2 
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Table 1 — Continued 



ID Galaxy RA (J2000) Dec Redshift L 4r a L radio h AOR Key c 



47 


UGC 5101 


09 35 51.7 


+61 21 11.3 


0.039 


11.96 


23.73 


4973056 


48 


3C 273 


12 29 06.7 


+02 03 08.6 


0.158 


12.83 


27.47 


4978176 


49 


Mrk 231 


12 56 14.2 


+56 52 25.2 


0.042 


12.51 


24.08 


4978688 


50 


Mrk 273 


13 44 42.1 


+55 53 12.7 


0.038 


12.09 


23.63 


4980224 


51 


Mrk 463 


13 56 02.9 


+18 22 19.1 


0.050 


11.80 


24.33 


4980736 


52 


Arp 220 


15 34 57.1 


+23 30 11.5 


0.018 


12.08 


23.34 


4983808 


53 


NGC 6240 


16 52 58.9 


+02 24 03.4 


0.024 


11.85 


23.92 


4985600 



"Infr ared luminosities are eithe r taken from Hurrah et all |2003) or calculated using the for- 
mula in lSanders fe Mirabell 1119961) . and converted to our cosmology. Units are the logarithm of 
the 1 — 1000/jm luminosity, in Solar luminosities (3.826xl0 26 W). 

b Obscrved-frame 1.4GHz luminosities, computed from the National Radio Astronomy Ob- 
servatory Very Large Array Sky Survey (NVSS) catalogs llCondon et al.lll998h , in units of log 
(W). The six objects without luminosities are not within the NVSS survey area. 

c Astronomical Observing Request Key number. 

d Redshift is tak en from the IR AS Point Source C atalog Redshift Survey (PSCz, 
ISaunders et al.ll2000h rather than from IWilman et al.l dl998l ) as our spectra are consistent with 
the PSCz redshift. 
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Table 2. Line ratio scaling factors for different extinction laws 



Line ratio Fluks a Li & Draine b Draine c Chiar & Tielens A d Chiar & Tielens B e 

[NcIII]A15.56 
[NcIIJA12.81 
[SIII]A18.71 



SIII 

SIV 


A:i:i.4K 
AH). 51 


SIII 

SIV' 


A18.71 
AH). 51 


SIII] 

OIV 


A33.48 
A25.89 
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NcV 


A12.81 
A14.32 


Noll 

'NcV 


A12.81 
A14.32 


NcV 


A24.32 



1.075 


1.052 


1.120 


0.881 


0.734 


0.684 


0.775 


0.634 




0.360 


0.498 


0.621 


0.506 


0.664 


0.373 


0.340 


0.478 


0.321 




0.134 


1.219 


1.149 


1.231 


1.148 


0.660 


1.163 


1.111 


1.202 


0.973 


1.000 


1.000 


0.990 


1.023 


0.968 


1.644 



Note. — Scaling factors assume an increase in the V band extinction of Ay = 30, in the 
form of a screen. 

HFluks et al] | | 1994ft 

iLi fc Drainej l|200lft 

torainj lll989ft 

jChiar fc Tielens! [120061) , ISM 

jChiar fc Tielens! j2006h . Galactic Center 



Table 3. Common emission lines 
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[bill] 
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ZLO lO 
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p. fox a 


97 ft 








91 fi 
Z1.0 


Q7 1 


11 .u 




9 1 ? 3 
zo.o 


Q7 1 
y ( . 1 


KA Q 




9^ ^ 
zo.o 


£ 9 


00188-0856 


<0.15 


0.75 


<0.26 


0.32 


4.67 


<0.18 


0.69 


0.73 


0.56:: 


<1.60 


<0.90 


1.39 






00397-1312 


0.35 


0.32: 


0.30 


<0.27 


4.41 


<0.20 


2.72: 


1.09:: 


1.80: 


<1.50 


<1.20 


<0.75 






01003-2238 




0.77 


0.21 


0.51: 


3.14 


<0.30 


1.31 


1.08 


<0.90 


<0.30 


<0.30 


<1.00 






03158+4227 


<0.60 


<1.40 


<1.20 


0.89 


5.78 


<1.10 


0.94:: 


2.86: 


1.71: 


<1.40 


<1.80 


<1.20 






03521+0028 


<0.42 


0.80 


<0.50 


0.70 


2.83 


<0.45 


1.27 


1.94 
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<0.48 


<0.90 


<0.51 






05189-2524 




3.22 


5.63 


1.67 


21.12 


17.53 


17.76 


3.15 


3.18 


11.73 


23.71 


<2.40 


<24.00 


11.85:: 


06035-7102 




2.05 


<0.39 


1.20 


6.98 


<0.48 


1.75: 


2.71 


3.00 


<0.81 


<3.00 


2.25 


4.61: 




06206-6315 




0.73 


0.21 


0.50: 


6.58 


2.30 


2.86 


1.36 


2.06 


2.03 


3.00 


<1.50 


3.53 




07598+6508 


<0.60 


<3.00 


<0.66 


0.78:: 


3.92 


<0.75 


2.45 


1.79 


<1.90 


<3.00 


<1.80 


<3.20 






08311-2459 




2.41 


8.70: 


1.68: 


24.29 


12.34 


22.53 


2.80 


11.43 


9.81 


26.30 


<4.40 


18.40: 




08572+3915 




0.43: 


<0.50 


<0.63 


7.18 


<0.75 


1.99: 


0.90 
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<5.40 


<6.00 


<2.50 


<12.00 


<25.00 


09022-3615 




6.32 


5.05 


6.02 


56.58 


<2.40 


40.00 


10.81 


24.73 


<3.60 


6.72 


11.16 


41.33 


95.83 


10378+1109 
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1.10 


<0.41 
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3.51 
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0.55:: 
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<2.40 
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10565+2448 




4.05 


<0.27 


2.53 
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6.42 
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<1.20 


1.85:: 


23.92 


43.40 


11095-0238 




1.80 
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1.19 


6.08 


<0.48 


1.89 
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1.22: 


<1.80 


<0.90 
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11119+3257 


<0.30 


0.52 


<1.50 


0.47:: 


2.97 


0.75 


1.98 


<2.60 


<1.80 


<0.90 


1.88:: 


<1.50 






12018+1941 


<0.20 


0.55: 


<0.19 


0.60: 


3.00 


<0.80 


0.35 


<1.80 


<0.70 


<0.10 


<0.63 


2.67: 






12071-0444 


0.61 


1.36 


1.46 


0.85 


5.25 


2.88 


5.09 


2.12 


1.86: 


3.71 


6.55 


<2.00 






12514+1027 


0.68 


0.27:: 


0.91 


0.50 


2.35 


1.94 


2.68: 


1.14:: 


1.21:: 


1.67: 


2.71: 








13120-5453 




6.89 


0.50: 


5.50 


150.04 


1.71 


18.46 


9.49 


19.18 


<20.00 


6.42: 


<5.00 


60.64 


107.13 


13218+0552 


<0.09 


0.44: 


<0.20 


0.25: 


1.07: 


<0.25 


<0.85 
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<3.40 
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0.40:: 
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0.35 
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Note. — Flux units are xlO 21 W cm 2 . A '-' indicates the line is outside the bandpass. Errors are of the order 10% or less, except for those fluxes 
marked with a ':' which are 20% and those marked with a '::' which are 30% and should be considered suspect. 
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Table 4. Unusual emission lines 
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Note. — Flux units arc xl0~ 21 W cm -2 . In this table, a '-' indicates that the line is not seen, 
however we do not quote upper limits due to the uncertain nature of the line IDs. Formal errors are 
all of the order 10%, except for ':' which are 20%, and '::' which arc 30%. 
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Fig. 1. — Short-High spectra of our sample. The numbers on the left hand side correspond 
to the ID numbers in column one of Table HJ Wavelengths have been shifted to the (optical) 
rest-frame . Some objects also show absorption at 13.7/im and 14.0/xm, caused by C 2 H 2 and 
HCN gas JLahuis et al.ll2007h . 
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2. — Short-High spectra (continued) 
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Fig. 3. — Short-High spectra (continued) 
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Fig. 6. — Long-High spectra of our sample. Note that the LH spectra are not arranged 
in order of ID number from Table [3] due to the wide dispersion in their continuum slopes. 
Wavelengths have been shifted to the (optical) rest-frame. Some objects also show an OH 
absorption feature at rest-frame 34.6//m. 
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Fig. 7. — Long-High spectra (continued) 
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Fig. 8. — Long-High spectra (continued) 
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Fig. 9. — Long-High spectra (continued). Note that the y axis scaling in this plot differs 
from that in Figures El [7] and El 
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Fig. 10. — Long-High spectra (continued). Note that the y axis scaling in this plot differs 
from that in Figures El [8] and M 
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Fig. 11. — Excitation diagnostic. The top panel identifies the ULIRGs using the numbers 
in column 1 of Table dj while the bottom panel shows the ULIRGs with error bars and 
3cr limits, as well as ancilliary data. ULIRGs with detections on both axes are plotted in 
black, whil e ULIRGs with lim its on one or both axes are p lotted in grey. Gre en symbols are 
starbursts (IVerma et al.ll2003l ) and red symbols are AGN (ISturm et al.l 120021 ). both with IR 



luminosities approximately in the range 10 10 - 1O 1L5 L . The dotted line is the fit to star- 
forming gala xies, while the da shed line is a fit to Seyfert galaxies, both with IR luminosities 
of < 10 10 L Q (IDale et al.l 120061 ) . The arrow indicates the effect on a points position if the V 
band extinction is increased by Ay = 30 (but see Table [2]). 
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Fig. 12. — [NeIII]A15.55/[NeII]A12.81 vs infrared luminosity. The small arrow indicates the 
effect on a points position if the V band extinction is increased by A v = 30. The horizontal 
bar on the bottom left indicates a 20% error on the IR luminosity. 
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Fig. 13. — [NeIII]A15.55/[NeII]A12.81 vs 1.4GHz luminosity. The small arrow indicates the 
effect on a points position if the V band extinction is increased by A v = 30. The horizontal 
bar on the bottom left indicates a 10% error on the 1.4GHz luminosity. 
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Fig. 14. — Total luminosity o f the [NeIHjA15.55 and [NeII]A12.81 lines vs IR luminosity. 
The solid line is the relation in lHo fc Ketd (120071 ). with the dashed lines indicating their lcr 
errors. ULIRGs with detections on both axes are plotted in black, while ULIRGs with limits 
on one or both axes are plotted in grey. The arrow indicates the effect on a points position 
if the V band extinction is decreased by Ay = 30. The horizontal bar on the bottom left 
indicates a 20% error on the IR luminosity. 
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Fig. 15.— The total luminosity of the [NeIII]A15.55 and [NeII]A12.81 lines vs the total 
luminosity of the PAH 6.2/im and 11.2/im features. ULIRGs with detections on both axes 
are plotted in black, while ULIRGs with limits on one or both axes are plotted in grey. The 
solid line indicates the best fit to the data, given in Equation [31 
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Fig. 16. — [NeV]A14.32/[NeII]A12.81 vs infrared luminosity. T he percentages on the y axis 
indicate the predicted AGN contribution to the IR luminosity (ISturm et al.ll2002l ). ULIRGs 
with detections on both axes are plotted in black, while ULIRGs with limits on one or both 
axes are plotted in grey. The small arrow indicates the effect on a points position if the V 
band extinction towards the line emitting regions is increased by Ay = 30. The horizontal 
bar on the bottom left indicates a 20% error on the IR luminosity. 
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Fig. 17. — [OIV]A25.89/[NeII]A12.81 vs infrared luminosity. Other data and symbols are the 
same as in Figure 
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Fig. 18.— [NeV]A25.89/[NeII]A12.81 vs IRAS 25/60 color. The green and red points are 
com parison starbu r sts an d A GN respectively, but un like Figure [H] these samples are taken 
from lBrandl et all J2006h and lWeedman et al.l (120051 ) . respectively. Other data and symbols 
are the same as in Figure [161 
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Fig. 19.— [NeV]A14.32/[NeII]A12.81 vs the equivalent width of the PAH 6.2/xm feature. The 
percentag es on the x axes in dicate the predicted contribution to the IR luminosity from an 
starburst (jArmus et al.l 120071 ). Other data and symbols are the same as in Figure fT6l 
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Fig. 20.— [OIV]A25.89/[NeII]A12.81 vs the equivalent width of the PAH 6.2/xm feature. 
Other data and symbols are the same as in Figure [191 
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Fig. 21.— [OIV]A25.89/[NeII]A12.81 vs the strength of the 9.7/im silicate feature. Other 
data and symbols are the same as in Figure [191 
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Fig. 22.— The total luminosity of the [NeIII]A15.55 and [NeII]A12.81 lines vs the total 
luminosity of the PAH 6.2/im and 11.2/zm features. The points are color coded according to 
the strength of the 9.7/im silicate feature. 
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Fig. 23.— [OIV]A25.89/[NeII]A12.81 vs the equivalent width of the PAH 6.2/im feature, with 
the points coded according to the strength of the 9.7 fiva. silicate feature. The color coding 
is the same as in Figure [22j Here we see that the lightly obscured and heavily obscured 
systems lie mostly on the left hand side of the plot, whereas moderately obscured systems 
lie mainly toward the right hand side of the plot. 



